Phylogenomics of Neogastropoda: the backbone hidden in the bush.

The molluscan order Neogastropoda encompasses over 15,000 almost exclusively marine species playing important roles in benthic communities and in the economies of coastal countries. Neogastropoda underwent intensive cladogenesis in early stages of diversification, generating a 'bush' at the base of their evolutionary tree, that has been hard to resolve even with high throughput molecular data. In the present study to resolve the bush, we use a variety of phylogenetic inference methods and a comprehensive exon capture dataset of 1,817 loci (79.6% data occupancy) comprising 112 taxa of 48 out of 60 Neogastropoda families. Our results show consistent topologies and high support in all analyses at (super)family level, supporting monophyly of Muricoidea, Mitroidea, Conoidea, and, with some reservations, Olivoidea and Buccinoidea. Volutoidea and Turbinelloidea as currently circumscribed are clearly paraphyletic. Despite our analyses consistently resolving most backbone nodes, three prove problematic: First, uncertain placement of Cancellariidae, as the sister group to either a Ficoidea-Tonnoidea clade, or to the rest of Neogastropoda, leaves monophyly of Neogastropoda unresolved. Second, relationships are contradictory at the base of the major 'core Neogastropoda' grouping. Third, coalescence-based analyses reject monophyly of the Buccinoidea in relation to Vasidae. We analysed phylogenetic signal of targeted loci in relation to potential biases, and we propose most probable resolutions in the latter two recalcitrant nodes. The uncertain placement of Cancellariidae may be explained by orthology violations due to differential paralog loss shortly after the whole genome duplication, which should be resolved with a curated set of longer loci.

Coupling DNA barcodes and exon-capture to resolve the phylogeny of Turridae (Gastropoda, Conoidea).

Taxon sampling in most phylogenomic studies is often based on known taxa and/or morphospecies, thus ignoring undescribed diversity and/or cryptic lineages. The family Turridae is a group of venomous snails within the hyperdiverse superfamily Conoidea that includes many undescribed and cryptic species. Therefore 'traditional' taxon sampling could constitute a strong risk of undersampling or oversampling Turridae lineages. To minimize potential biases, we establish a robust sampling strategy, from species delimitation to phylogenomics. More than 3,000 cox-1 "barcode" sequences were used to propose 201 primary species hypotheses, nearly half of them corresponding to species potentially new to science, including several cryptic species. A 110-taxa exon-capture tree, including species representatives of the diversity uncovered with the cox-1 dataset, was build using up to 4,178 loci. Our results show the polyphyly of the genus Gemmula, that is split into up to 10 separate lineages, of which half would not have been detected if the sampling strategy was based only on described species. Our results strongly suggest that the use of blind, exploratory and intensive barcode sampling is necessary to avoid sampling biases in phylogenomic studies.

Systematics and Distributions of Upper Bathyal Species in Bathyancistrolepis, a Deep-Sea Whelk Genus Endemic to the Northwest Pacific (Gastropoda: Buccinidae).

The deep-sea buccinid snail genus Bathyancistrolepis is redefined based on the reconstruction of a molecular phylogeny and morphological examination of shell and radular characters. This genus is distinguished from other genera of the subfamily Parancistrolepidinae with a combination of shell traits, including (1) a low spire, (2) sharp, carinate spiral cords or keels and (3) a long, curved siphonal canal, but not with a difference in radular morphology as suggested by previous authors. Three allopatric or parapatric species are recognized in the upper bathyal (447-2057 m) waters around Japan and Taiwan: B. tokoyodaensis from off Hokkaido to Sagami Bay in the Northwest Pacific, B. trochoidea off Kumano-nada to Miyazaki in the Northwest Pacific and along Nansei Islands in the East China Sea, and B. taiwanensis sp. nov. in the South China Sea. These species bear large paucispiral protoconchs that are indicative of benthic early development without a pelagic larval period, and hence low dispersal capability. Seafloor topography seems to have acted as a barrier for their dispersal; the range of B. tokoyodaensis supports the previous finding that Izu Peninsula delimits westward distribution of bathyal gastropod species of boreal origins.

Neogastropod (Mollusca, Gastropoda) phylogeny: A step forward with mitogenomes.

The Neogastropoda (Mollusca, Gastropoda) encompass more than 15,000 described species of marine predators, including several model organisms in toxinology, embryology and physiology. However, their phylogenetic relationships remain mostly unresolved and their classification unstable. We took advantage of the many mitogenomes published in GenBank to produce a new molecular phylogeny of the neogastropods. We completed the taxon sampling by using an in-house bioinformatic pipeline to retrieve mitochondrial genes from 13 transcriptomes, corresponding to five families not represented in GenBank, for a final dataset of 113 taxa. Because mitogenomic data are prone to reconstruction artefacts, eight different evolutionary models were applied to reconstruct phylogenetic trees with IQTREE, RAxML and MrBayes. If the over-parametrization of some models produced trees with aberrant internal long branches, the global topology of the trees remained stable over models and softwares, and several relationships were revealed or found supported here for the first time. However, even if our dataset encompasses 60% of the valid families of neogastropods, some key taxa are missing and should be added in the future before proposing a revision of the classification of the neogastropods. Our study also demonstrates that even complex models struggle to satisfactorily handle the evolutionary history of mitogenomes, still leading to long-branch attractions in phylogenetic trees. Other approaches, such as reduced-genome strategies, must be envisaged to fully resolve the neogastropod phylogeny.

Venom Diversity and Evolution in the Most Divergent Cone Snail Genus Profundiconus.

Profundiconus is the most divergent cone snail genus and its unique phylogenetic position, sister to the rest of the family Conidae, makes it a key taxon for examining venom evolution and diversity. Venom gland and foot transcriptomes of Profundiconus cf. vaubani and Profundiconusneocaledonicus were de novo assembled, annotated, and analyzed for differential expression. One hundred and thirty-seven venom components were identified from P. cf. vaubani and 82 from P. neocaledonicus, with only four shared by both species. The majority of the transcript diversity was composed of putative peptides, including conotoxins, profunditoxins, turripeptides, insulin, and prohormone-4. However, there were also a significant percentage of other putative venom components such as chymotrypsin and L-rhamnose-binding lectin. The large majority of conotoxins appeared to be from new gene superfamilies, three of which are highly different from previously reported venom peptide toxins. Their low conotoxin diversity and the type of insulin found suggested that these species, for which no ecological information are available, have a worm or molluscan diet associated with a narrow dietary breadth. Our results indicate that Profundiconus venom is highly distinct from that of other cone snails, and therefore important for examining venom evolution in the Conidae family.

Exon-Capture-Based Phylogeny and Diversification of the Venomous Gastropods (Neogastropoda, Conoidea).

Transcriptome-based exon capture methods provide an approach to recover several hundred markers from genomic DNA, allowing for robust phylogenetic estimation at deep timescales. We applied this method to a highly diverse group of venomous marine snails, Conoidea, for which published phylogenetic trees remain mostly unresolved for the deeper nodes. We targeted 850 protein coding genes (678,322 bp) in ca. 120 samples, spanning all (except one) known families of Conoidea and a broad selection of non-Conoidea neogastropods. The capture was successful for most samples, although capture efficiency decreased when DNA libraries were of insufficient quality and/or quantity (dried samples or low starting DNA concentration) and when targeting the most divergent lineages. An average of 75.4% of proteins was recovered, and the resulting tree, reconstructed using both supermatrix (IQ-tree) and supertree (Astral-II, combined with the Weighted Statistical Binning method) approaches, are almost fully supported. A reconstructed fossil-calibrated tree dates the origin of Conoidea to the Lower Cretaceous. We provide descriptions for two new families. The phylogeny revealed in this study provides a robust framework to reinterpret changes in Conoidea anatomy through time. Finally, we used the phylogeny to test the impact of the venom gland and radular type on diversification rates. Our analyses revealed that repeated losses of the venom gland had no effect on diversification rates, while families with a breadth of radula types showed increases in diversification rates, thus suggesting that trophic ecology may have an impact on the evolution of Conoidea.

Quid est Clea helena? Evidence for a previously unrecognized radiation of assassin snails (Gastropoda: Buccinoidea: Nassariidae).

The genus Clea from SE Asia is from one of only two unrelated families among the megadiverse predatory marine Neogastropoda to have successfully conquered continental waters. While little is known about their anatomy, life history and ecology, interest has grown exponentially in recent years owing to their increasing popularity as aquarium pets. However, the systematic affinities of the genus and the validity of the included species have not been robustly explored. Differences in shell, operculum and radula characters support separation of Clea as presently defined into two distinct genera: Clea, for the type species Clea nigricans and its allies, and Anentome for Clea helena and allies. A five-gene mitochondrial (COI, 16S, 12S) and nuclear (H3, 28S) gene dataset confirms the placement of Anentome as a somewhat isolated offshoot of the family Nassariidae and sister to the estuarine Nassodonta. Anatomical data corroborate this grouping and, in conjunction with their phylogenetic placement, support their recognition as a new subfamily, the Anentominae. The assassin snail Anentome helena, a popular import through the aquarium trade so named for their voracious appetite for other snails, is found to comprise a complex of at least four species. None of these likely represents true Anentome helena described from Java, including a specimen purchased through the aquarium trade under this name in the US and one that was recently found introduced in Singapore, both of which were supported as conspecific with a species from Thailand. The introduction of Anentome "helena" through the aquarium trade constitutes a significant threat to native aquatic snail faunas which are often already highly imperiled. Comprehensive systematic revision of this previously unrecognized species complex is urgently needed to facilitate communication and manage this emerging threat.

A multilocus molecular phylogeny of Fasciolariidae (Neogastropoda: Buccinoidea).

The neogastropod family Fasciolariidae Gray, 1853 - tulips, horse-conchs, spindles, etc., comprises important representatives of tropical and subtropical molluscan assemblages, with over 500 species in the subfamilies Fasciolariinae Gray, 1853, Fusininae Wrigley, 1927 and Peristerniinae Tryon, 1880. Fasciolariids have had a rather complicated taxonomical history, with several genus names for a long time used as waste baskets to group many unrelated species; based on shell characters, recent taxonomic revisions have, however, began to set some order in its taxonomy. The present work is the first molecular approach to the phylogeny of Fasciolariidae based on a multigene dataset, which provides support for fasciolariids, an old group with a fossil record dating back to the Cretaceous. Molecular markers used were the mitochondrial genes 16S rRNA and cytochrome c oxidase subunit I, and the nuclear genes 18S rRNA, 28S rRNA and histone H3, sequenced for up to 116 ingroup taxa and 17 outgroups. Phylogenetic analyses revealed monophyly of Dolicholatirus Bellardi, 1884 and Teralatirus Coomans, 1965, however it was not possible to discern if the group is the sister clade to the remaining fasciolariids; the latter, on the other hand, proved monophyletic and contained highly supported groups. A first split grouped fusinines and Pseudolatirus Bellardi, 1884; a second split grouped the peristerniine genera Peristernia Mörch, 1852 and Fusolatirus Kuroda and Habe, 1971, while the last group comprised fasciolariines and the remaining peristerniines. None of these clades correspond to the present-day accepted circumscription of the three recognized subfamilies.

THE INDO-PACIFIC GEMMULA SPECIES IN THE SUBFAMILY TURRINAE: ASPECTS OF FIELD DISTRIBUTION, MOLECULAR PHYLOGENY, RADULAR ANATOMY AND FEEDING ECOLOGY.

The biology, feeding ecology and phylogenetic relationships of marine snails in the family Turridae remain poorly understood. Here we report our study on four deep-water species in the genus Gemmula, a major group in this family. The four species G. speciosa (Reeve 1843), G. sogodensis (Olivera 2005), G. kieneri (Doumet 1940) and G. diomedea (Powell 1964) were collected at five different sites in the Philippines, and their pattern of distribution in the sites, their feeding behaviour as well as their phylogenetic relationships with each other and with other members of the subfamily Turrinae were investigated. The radular morphology (of two Gemmula species) and potential prey (for one Gemmula species) were also examined. Actual feeding observations were also conducted for Gemmula speciosa and compared with two turrids from other genera.All four Gemmula species showed strikingly different patterns of distribution; each species was found to be relatively much more abundant at one site but not at the other sites. Molecular phylogenetic analysis based on 16S sequences correlated with previously reported 12S sequences and revealed that the four species all belong to a well-supported Gemmula clade within the subfamily Turrinae; and that this clade appeared more closely related to the clades Xenuroturris, Turris and Lophiotoma than to the other clades in the subfamily (i.e., Turridrupa, Unedogemmula and Polystira). Morphological analysis of the radula of both G. speciosa and G. sogodensis revealed that the radulae of the two species were similar but differed from the other turrids, Lophiotoma acuta and Unedogemmula bisaya, by the absence of central teeth, consistent with the separation of the Gemmula clade from the Lophiotoma and Unedogemmula clade.To identify the polychaete group that is targeted as prey by species of Gemmula, analysis of regurgitated food fragments was made; phylogenetic analysis of an mtCOI gene fragment that was PCR-amplified from the regurgitated tissue of one specimen (G. diomedea) indicated close affinity of the prey to the terebellid polychaete Amphitritides. Specimens of Gemmula speciosa, when challenged with the terebellid polychaete Loimia sp., were observed to attack the worm suggesting that Gemmula species feed on terebellid polychaetes. Lophiotoma acuta were also observed to feed on the same species of terebellid but were usually group-feeding in contrast to the solitary feeding of G. speciosa. Unedogemmula bisaya did not feed on the terebellid which also supports the separation of the Gemmula and Unedogemmula clade.Two lines of proof (i.e. the molecular phylogenetic analysis and the feeding challenge) supporting the toxin homology findings previously reported, provide consistent evidence that Gemmula is a distinct clade of worm-hunting Turrinae that feeds on Terebellidae.

Accessing novel conoidean venoms: Biodiverse lumun-lumun marine communities, an untapped biological and toxinological resource.

Cone snail venoms have yielded pharmacologically active natural products of exceptional scientific interest. However, cone snails are a small minority of venomous molluscan biodiversity, the vast majority being tiny venomous morphospecies in the family Turridae. A novel method called lumun-lumun opens access to these micromolluscs and their venoms. Old fishing nets are anchored to the sea bottom for a period of 1-6months and marine biotas rich in small molluscs are established. In a single lumun-lumun community, we found a remarkable gastropod biodiversity (155 morphospecies). Venomous predators belonging to the superfamily Conoidea (36 morphospecies) were the largest group, the majority being micromolluscs in the family Turridae. We carried out an initial analysis of the most abundant of the turrid morphospecies recovered, Clathurella (Lienardia) cincta (Dunker, 1871). In contrast to all cDNA clones characterized from cone snail venom ducts, one of the C. cincta clones identified encoded two different peptide precursors presumably translated from a single mRNA. The prospect of easily accessing so many different morphospecies of venomous marine snails raises intriguing toxinological possibilities: the 36 conoidean morphospecies in this one net alone have the potential to yield thousands of novel pharmacologically active compounds.

Alpha-conopeptides specifically expressed in the salivary gland of Conus pulicarius.

To date, studies conducted on cone snail venoms have attributed the origins of this complex mixture of neuroactive peptides entirely to gene expression by the secretory cells lining the lumen of the venom duct. However, specialized tissues such as the salivary glands also secrete their contents into the anterior gut and could potentially contribute some venom components injected into target animals; evidence supporting this possibility is reported here. Sequence analysis of a cDNA library created from a salivary gland of Conus pulicarius revealed the expression of two transcripts whose predicted gene products, after post-translational processing, strikingly resemble mature conopeptides belonging to the alpha-conotoxin family. These two transcripts, like alpha-conotoxin transcripts, putatively encode mature peptides containing the conserved A-superfamily cysteine pattern (CC-C-C) but the highly conserved A-superfamily signal sequences were not present. Analysis of A-superfamily members expressed in the venom duct of the same C. pulicarius specimens revealed three putative alpha-conotoxin sequences; the salivary gland transcripts were not found in the venom duct cDNA library, suggesting that these alpha-conotoxins are salivary gland specific. Therefore, expression of conotoxin-like gene products by the salivary gland could potentially add to the complexity of Conus venoms.

Morphological proxies for taxonomic decision in turrids (mollusca, neogastropoda): a test of the value of shell and radula characters using molecular data.

The state of the art of turrid (=Turridae s. l.) systematics is that shells - when they include the protoconch - are reliable species-level identifiers, but inadequate proxies for allocation to genera or subfamilies. Generally, the radula is used for allocation to a (sub)family, but the hypothesis that the radula is a more adequate proxy than the shell for relationships has not yet been tested by molecular data. Species of Xenuroturris may have drastically different radulae, with either "semi-enrolled" or "duplex" marginal teeth, although their shells are very similar or even almost indistinguishable. Molecular data confirm that specimens with different types of radulae constitute different species, but two species of a pair with respectively semi-enrolled and duplex teeth end up being not closely related. However, it is still unresolved whether species with semi-enrolled (= Iotyrris ) and duplex teeth (= Xenuroturris ) form two supported monophyletic clades. Iotyrris devoizei n.sp. and I. musivum n.sp. are described from Vanuatu, where they occur sympatrically with I. cingulifera and Xenuroturris legitima .